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PREFACE 

Invariant  imbedding  techniques  have  proven  useful  in  the  compu¬ 
tational  solution  of  problems  in  neutron  transport  and  radiative  trans¬ 
fer.  This  Memorandum  discusses  an  application  of  these  techniques  to 
a  problem  in  rarefied  gas  dynamics  and  should  be  of  interest  to  aero- 
dynamicists,  physicists,  and  applied  mathanaticians. 
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SUMMARY 


This  Memorandum  re-examines  the  application  of  invariance  tech¬ 
niques  to  plane  Couette-flow  problems  and  presents  several  different 
approaches  to  the  derivation  of  the  invariant- imbedding  equations. 

It  is  shown  how  both  linear  and  nonlinear  collision  operators 
may  be  considered,  for  both  discrete  and  continuous  versions  of  the 
original  Boltzmann  formulation. 
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SYMB0L3 

velocity  vector,  defined  by  three  velocity  components 

Ui  ,  Vi  ,  Wi^ 

U  V  w 

distribution  function  for  particles  with  velocity 
vector  (?i^ 

th 

collision  operator  for  particles  In  1  velocity 
state 

linearized  collision  operator 

linearized  collision  operator  acting  on  i 

linearized  collision  operator  acting  on  i 

reflection  function 

nonlinear  reflection  function 

defined  by  Eq.  (20) 

one>dlmenslonal  reflection  function 

molecular  velocities  for  particles  In 

velocity  components  In  discrete  velocity  space 

molecular  velocities  for  particles  In  i 
distance  along  normal  to  a  moving  plate 
spacing  between  plates 

source  function  for  molecules  In  1^^  state  entering 
at  X  "  z 

perturbation  distribution  function 

perturbation  distribution  function  for  particles  with 
velocity  components  In  the  direction  of  Increasing  : 

perturbation  distribution  function  for  particles  with 
velocity  components  In  the  direction  of  decreasing  x 


. . . . . . . . 
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I,  INTROroCTION 

In  an  earlier  paper, che  concepts  of  invariant  imbedding  were 
applied  to  a  simple  problem  in  rarefied  gas  dyttamics,  that  of  steady, 
linearized  shear  flow  between  two  infinite  plane  walls  under  the  as¬ 
sumptions  of  Krook  scattering  and  diffuse  reflection  from  both  walls. 
Particle-counting  techniques  were  applied  directly  to  the  Krook  col¬ 
lision  operator,  reflection  functions  were  introduced,  and  the  charac¬ 
teristically  nonlinear  integrc-dif ferential  equation  of  the  initial- 

(2) 

value  tJTpe  ^or  the  reflection  function  was  derived.  Subsequent  work 
has  concentrated  on  the  computational  solution  of  that  integral  equa¬ 
tion  as  well  as  on  the  unsteady  version  of  this  same  linearized  shear 
flow  and  on  the  effect  of  specular-diffuse  molecular  reflection  from 
boundaries . 

In  this  Memorandum,  we  briefly  re-examine  the  application  of  in¬ 
variance  techniques  to  plane  Couette-flow  problems  and  present  several 
alternate  prescriptions  for  the  derivation  of  the  invariant- imbedding 
equations  di.rectly  from  the  Boltzmann  formulation*  Our  objective  here 
is  not  to  present  any  new  results  but  rather  to  review  and  extend 
several  techniques  which  have  already  been  developed. 

It  is  clear  by  now  that  a  more  formal  and  mechanical  approach  than 
particle  counting  Is  desirable  for  the  derivation  of  invariance  equations. 
Particle  counting,  unless  performed  with  great  care,  can  lead  to  errors 
in  problems  involving  complex  geometries  and  nonlinear  scattering  models. 

A  rigorous  treatment  of  the  equations  of  invariant  imbedding  has  recent- 

(3) 

ly  been  given  by  Bailey,  who  considered  a  number  of  related  linear 
prrhlems  in  transport  theory.  In  particular,  ha  has  re-examined  the 
earlier  analyses  for  spherical  geometries  and  concludes  that  perturba¬ 
tion  techniques  are  the  most  direct  method  of  approach. 

The  invariant- imbedding  approach  involves  the  calculation  of  the 
effect  of  a  variation  of  the  spacing  between  the  two  boundaries  for  a 
fixed  Incoming  distribution  function  of  a  particle  at  one  boundary.  We 

shall  consider  two  different  approaches  to  the  problem;  the  first,  that 
(4) 

>f  Wing,  '  is  suitable  for  linear  problems,  while  the  second,  that  of 
Bellman  and  Kslaba,^^^  is  suitable  for  both  linear  and  nonlinear  problems. 
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II.  ANALYSES 


METHOD  I 

We  consider  a  linearized  Couette  flow  between  two  infinite  pla> 
nar  walls.  The  Boltzmann  equation  for  this  flow  becomes 


hi 

U  -  (U.  V.  w> 


U^) 


(1) 


The  notation  of  Ref.  1  is  employed,  where  u,  v,  and  w  are  dimension¬ 
less  molecular  velocities,  4  is  a  perturbation  distribution  function, 

X  is  distance  in  units  of  mean  ftcie  path,  and  L($)  is  a  linearized 
collision  operator.  It  was  shown  in  Ref.  2  how  it  is  possible  to  un¬ 
couple  the  effect  of  the  boundary- surface  interaction  by  dealing  with 
a  standard  s^  of  boundary  conditions  corresponding  to  one  wall  being 
a  sink  for  itolecules  and  the  other  wall  being  a  molecular  source.  This 
uncoupling  is  only  valid  for  linearized  flows.  The  solution  of  this 
"standard"  problem  can  then  be  used  to  construct  solutions  which  include 

the  effects  of  nondiffuse  surface-particle  interaction.  The  standard 

* 

boundary  conditions  are 


4'’'(u,  V,  w,  f(u  ,  V  ,  w  ),  0) 


o  o  o 


ooo  ooo  ooo 


(sink  at  X  •  0) 

(2) 

(arbitrary  sor.rce 
of  molecules  at 
X  ■  z) 


where  we  have  chosen  to  indicate  explicitly  the  dependence  of  the  4's 
on  f(u,  V,  w) .  We  seek  a  function  R(u,  v,  w,  u^,  w^) ,  called  a 

reflection  function,  such  that 


The  subscripted  varlsb’.rs  u  ,  v  ,  w  refer  to  the  velocity 
space  of  particle®  with  velocity  components  in  the  direction  of  de¬ 
creasing  X,  ''•r  cnoea  particles  in  i  ,  while  the  unsubccripted  vari¬ 
ables  V,  V  refer  to  the  velocity  space  of  particles  with  velocity 
cn:tpcnents  in  the  dlrectlor.  of  increasing  x,  or  those  particles  in 
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o  j,-*  «■+« 

i  (u,  V,  w,  f(u,  V,  w),  *)  «  R(u,  V,  W,  U  ,  V  ,  w  ,  z) 

v'  j  V  0  0  0 

•  C9  ^00 


£(u  ,  V  ,  w  )  du  dv  dw 
o  o  o  o  o  o 


(3) 


We  consider  the  same  problem  as  that  given  by  Eq.  (3)  and  boundary 
conditions  of  Eq.  (2)  on  the  interval  0  ^  x  ^  z  +  A.  Then 

$^(u,  V,  w,  f(u  ,  V  ,  w  ,  0)  «  0 

(4) 

*  (“o'  “o'  "o'  ((“o'  “o'  Xo).  '  +  *)  ■ 


and 


,0  A+®  «+* 


(u,  V,  w,  f(u^,  v^,  w^),  z  +  A)  -  J  J  J  R(u,  V,  w,  u^,  v^,  w^,  z  +  A) 


^m^co  -a 


f(u  ,  V  ,  w  )  du  dv  dw 

O  O  O  0  0  0 


(5) 


From  Eq.  (1),  we  can  write 
$'*'(z  +  A,  f)  * 

§''’(z,  f)  +  “  (4^(u,  V,  w,  f,  z),  v^,  w^,  f,  z)  +  o(A))  (6) 


where 


lim 

A-K) 


o(A} 

A 


0 


and  f  s  f(u  ,  v  ,  w  ).  We  now  observe  that  if  we  consider  this  p.ob- 
o  o  o 

lem  on  the  interval  0  ^  x  ^  z  we  may  write 


4- 


nr 

•90  *00  •ao 


i  (u,  V,  w,  f,  *)  - 


R(u,  V,  w,  u^,  v^,  w^,  z) 


f  (u  ,  V  ,  w  )  du  dv  dw 
o  o  o  o  o  o 


where  f  (u  ,  v  ,  w  )  is  the  source  function  for  the  new  problem. 
''  o  o 


f  (u  ,  V  ,  W  )  ■  $  (u  ,  V  ,  W  ,  fi  2) 
o  o  o  o  o  o 


f  (u  ,  V  ,  w  ,  f,  z  +  A)  -  —  L  (z  +  A;  +  o(A) 
o  o  o  u 

o 


f(u  ,  V  ,  w  )  -  -  (z  +  A)  +  o(A) 

o  o  o  u 

o 


Thus,  from  Eqs.  (5),  (6),  And  (8),  we  may  write 

rf  I  R(u,  V,  w,  u  ,  V  ,  w  ,  z  +  A)f(u  ,  V  ,  V  )  du  dv  dw 
**«  o  o  '  '  o’  o*  o  o  o  o 

■  j  j  r  ”•  “o'  “o'  “o' 


*ao  •«  •« 


“o'  "o> 


X  du  dv  dw  +  —  AL*'(z)  +  o(A) 
o  o  o  u 


We  note  that 


AL*(z  +  A)  ■  AL*(z)  +  o(A) 


Equation  (9)  contains  the  gist  of  the  invariant -imbedding  argument. 

The  Krook  equation  possesses  the  property  that  the  perturbation 
distribution  func  .ion  and  the  linearized  collision  operator  are  reduc¬ 
ible  to  a  dependence  on  the  single  velocity  variable,  u,  for  linearized 


If  the  linearized  collision  operator  L  is  specified,  then  an 
integro* differential  equation  for  R(u,  v,  w,  u  ,  v  ,  w  )  is  recovered. 

o  o  o 
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9te£dy  and  unataady  shear  flows  (boundary  conditions  permitting).  This, 
of  course,  Is  a  major  virtue  of  this  scattering  model.  Consequently, 
li'tegrations  over  v^  and  are  performed  explicitly,  and  a  significant 
reGuctiou  In  the  number  of  Independent  variables  Is  obtained.  An  equa¬ 
tion  Is  obtained  for  the  one- dimensional  reflection  function  R(u,  u^,  t) 
by  -'ochlnlng  Eqs.  (9)  and  (10). 


.+  ^1  L+  JL  r  If 

'  Llll.lon,  1 


,rf>  ”U 

1  e  °  S  du 


collisions 


r  •>  ^ 

1C®  2  ,  ^  -u 

“  "  1  *  du-*^je  *■ 

V. 


$  du 


Introducing  Eq.  (10)  Into  Eq.  (?),  we  obtain 


r  R(u,  u  ,  *  +  A)f(u  )  du  -  f  R(u,  u  ,  a)-"  f(u  )  +  4“ 
o  ooj  o  ou 


)  du 
o  o 


*  2  ^ 

+  "7=  “  A  r  e  I  Rlu',  u',  t)f(u*)  du*  du* 

yn  u  J_  J_  o’  '  '  o  o 


+  7=  -  A  I  f(u’)e  °  du* 
Vn  u  J  o  o 

▼  «fl0 


where  u  ,  u',  and  u*  are  dumy  variables. 

o  o 
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This  relationship  roust  hold  for  all  source  functions  f (u) .  We 
select  Che  particular  value  corresponding  to  a  Dirac-delta-functlon 

source ,  i ■ e  , 


f(u)  ■  o(u  -  u^) 


Therefore, 


R(u,  u^,  T  +  A) 


R(u,  u^,  z)  +  “  u^,  z) 

o 


r®  1 2 

•  Ww  f  ‘*''0 

-«  oi  ''  o 


-  4  f  R(U.  u;.  z)  jfe  °  dul 


t  “o- 


«  ,  2  -u^ 

rt 


Letting  A  -•  0,  and  Introducing 


1  1  ”% 

R(U.  u^,  t)  ‘  °  !>(u.  |u^|,  z) 


Ln  order  to  better  obtr/n  a  synmetric  ic’^u  yislds 


(12) 


iL 
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^  (u,|  u^l,  z)  +  -j— r  (u,  |u^|.  z)  +  -  (u,  |u^j,  z) 


,  -u 

1  I  e 


-u'" 

,  ,®  o 
1  I  e 


+  --  1  - S(u,  u'  ,  z)  du' 

ttJ  iu'i  '  o  o 


o  o 


,2  „ 
•“I  _,..2 


M  r  ‘iu'i!  u' i“o!’  i“oi‘  “*“0 


This  equation  was  derived  in  Ref.  1  by  particle-countitig  tech¬ 
niques  . 

METHOD  II 

A  direct  and  powerful  approach  to  the  development  cf  the  invariant- 
imbedding  equations  from  the  Boltzmann  formulation  uses  simple  pertur¬ 
bation  techniques  as  the  principal  cool.  We  apply  it  here  to  a  discrete 
•'del  of  the  Boltzmann  equation  and  then  indicate  the  results  for  con¬ 
tinuous  models. 

Dlacrete  Model 

We  introduce  a  discrete  model  of  the  Boltzmann  equation: 


c  “  ^2  *“  ^l’  ^2  **’ 


dx  “  "*^i^^l’  ^2’  *“  ^N’  ^1’  ^2  ***  V 
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Tlie  i^’s  are  related  to  discrete  versions  of  the  collision  opera¬ 
tor  and  are  not  restricted  solely  to  linear  operators  on  the  f^’s. 
Discrete  models  of  transport  equations  have  been  considered  for  many 
years  in  neutron  and  radiative  transfer,  but  it  is  only  in  recent 
years  that  attempts  have  been  made  to  utilize  them  in  kinetic  theory. 
We  note  that  the  subscript  i  refers  to  a  particular  value  of  the 
velocity  vector,  C^,  defined  by  =  ^'^iu*  '^ly'  ^i 

to  the  distribution  function  of  molecules  having  positive  u  compo¬ 
nents,  and  the  refer  to  molecules  having  negative  u  components 

where  u  is  positive  in  the  direction  of  increasing  x. 

The  discrete  model  arises  from  the  consideration  of  a  discrete 
set  of  molecular  velocities,  rather  than  a  continuous  set.  Finite 
volumes  in  velocity  space  *•  j  utilized,  rather 

than  infinitesmal  elements  dc*  =  du  dv  dw;  integrals  are  replaced  by 
sums,  and  a  distribution  function  corresponding  to  each  of  the  states 
is  sought . 

We  consider  boundary  conditions  resulting  from  an  arbitrary 
source  of  molecules  at  x  “  z  and  a  sink  for  molecules  at  x  ■»  0.  A 
somewhat  improbable  physical  interpretation  of  these  boundary  condi¬ 
tions  is  that  the  wall  at  x  ■  z  evaporates  molecules  into  the  flow 
and  possesses  a  unit  sticking  coefficient,  while  the  wall  at  x  ■>  0 
condenses  all  molecules  which  impinge  on  it.  For  linear  problems, 
as  discussed  earlier,  this  standard  set  of  boundary  condltlo.<.a  is 
sufficient  to  account  for  other  more  interesting  surface-particle 
interactions. 

These  boundary  conditions  are 


(15) 


Let  the  solution  of  Eq.  (14)  on  the  interval  0,  z,  subject  to  the 
source  at  x  ■  z,  be  indicated  by  f^,  f^^. 

Similarly,  let  the  solution  of  Eq.  (15)  on  the  Interval  0,  z  f  A, 
subject  to  the  source  at  x  ■  z  +  A,  be  Indicated  by  F^,  f”. 
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Thus  we  ttsy  write 


'**1  +  +  + 


.  ^  f" 

N’  1’  2 


dx  ’  ^2  “•  ^N*  ^1*  ^2  *"  V 


where  1  runs  from  I  Co  N  with  boundary  conditions 


£j(a)  -  0 


£'(«)  -  y; 


'“l  +  -  -  . 

dT-V^-Fj'- V  Fp  r^.-.F^) 


jji-  -K-(F+  F*  ...  fJ.  F;.  F-  ...  F-) 


where  the  is  the  same  collision  operator  i.4  the  k*,  where 

F^(0)  -  0 


f‘(*  +  a)  -  Yi 


We  seek  perturbation  solutions  to  Eq.  (18)  in  the  form 


F^  -  f  +  A  •  r^  +  o(A) 


-  fj^  +  A  •  Wj^  +  o(A) 


(20) 


The  boundery  condltione  on  r^,  can  be  obtained  as  follows: 


r^(0)  -  0  at  X  "  0 


(21a) 


At  X  «  z>  we  expand  in  a  Taylor  series  to  obtain 


f'(z)  -  F‘(z  +  A)  -  a  •  ^  (e  +  A)  +  0(A) 


(21b) 


i~A^)  +  A  •  wj(x)  -  yI  -  A  •  ^ 


(a)  +  o(A) 


w"(x)  -  k'(f+(z),  f*iz)  ...  fj(£),  f"(£),  f“(E)  ...  f‘(E))  (21c) 


The  linear  differential  equations  for  the  perturbations  r^»  v~  are 
obtained  by  subswitution  of  Eq.  (20)  into  Eq.  (18). 


dr,  dk. 


ok  . 


i  V  i  +  j.  V  i  - 

ir  -  Z  +  Z,  "j 

j-i  j-i 


.1  8k 


i  V  i  +  V  i 

""  -  Z  ■'Z 


j-1  j 


The  solution  to  this  set  of  equations,  subject  to  boundary  con¬ 
ditions  in  Eqs.  (21),  is 
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N 


df 


*(x)  .  ^  —i  .  k'(f^<»).  fj(,)  ... 


N 


f^C*)  ...  f‘(*)  ...  f'(*)) 


(23) 


r 

■  1  tjc*)  ••■  f;( 


J-1 


»)) 


which  can  be  verified  by  aubstitutlon. 

At  X  ■  *  +  A,  a  Taylor  aeries  expansion  results  In 

+  +  AaF+Cs) 

F^(*  +  A)  -  F^(*)  +  — - +  o(A) 


+  Arl’Cs) 


+  Ak^(f{\s),  £^(s) 


^u(*)»  ■*■  °(^) 


(24) 


At  this  point  in  the  analysis  it  is  appropriate  to  introduce  a 
reflection  function  that  must  differ  essentially  from  the  reflection 
function  of  the  linear  problem  considered  earlier.  In  the  linear 
problem,  it  was  assumed  that  the  distribution  function  for  particles 
lmplx.glng  on  the  boundary  at  x  -  z  was  a  linear  functional  of  the  dis¬ 
tribution  function  for  particles  emitted  from  that  boundary.  The 
definition  of  the  reflection  function  in  Eqs.  (3)  and  all  the  succeed¬ 
ing  analysis  was  engendered  by  the  linearity  of  the  collision  operator. 
Thus  ..t  was  possible  to  consider  a  delta- function  incident  flux  of 
molecules  and  construct  from  this  the  solution  for  any  Incident  flux. 

In  this  nonlinear  vision  of  the  same  pioblem  it  is  no  longer  possible 
to  uncouple  the  ei.  of  a  particular  source  distribution  at  the 
boundary  at  x  -  z  from  the  resulting  flux  of  particles  incident  upon 
that  boundary.  We  introduce  the  function  R*(z,  vj  • • •  Y^)  which 
is  the  distribution  function  of  particles  incident  upon  the  boundary 
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aL  X  -  z  due  to  the  particular  distribution  of  incident  particles  at 

X  =  2  emitted  by  that  boundary,  V-i  ♦  Yn  Y«' 

1  i  N 

Thus 


Yi . 


r  ^2  • 


y^)  £^(2.) 


Rj_(z  +  A,  Y2 


vjj)  -  F^(2  +  A) 


and 


lim 

A-O 


F^Cz  +  A)  -  f^(2) 


The  resulting 


equation  for  R^Cz)  is  then 


dR^(x,  Y2  • ' '  Yjj) 


X  k^(R^(z),  R^Cz)  ...  ys),  Yj,  Y2  •••  Y^) 

+k^(Rj(z),  R^(z)  ...  R^(z),  vj,  Y2  •••  Y^) 

(25) 


Continuous  Model 


If  a  continuous  version  of  Eqs.  (16)  and  (17)  is  considered, 
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^  (c‘,  x)  -  -k'cc"^) 


£"^(0^,  0)  *  0 


(26) 


f‘(C  ,  2)  -  y(C  ) 

“*±  + 
where  the  K(C  )  are  functionals  of  the  distribution  functions  r  ,  and 

C  takes  on  all  values  of  the  velocity  vector  possessing  a  component 

"H- 

in  the  direction  of  decreasing  x,  while  C  takes  on  all  values  posses¬ 
sing  a  component  in  the  direction  of  increasing  x.  The  results  of  the 
previous  section  taay  be  extended  formally  to  give 


(y(c),  c,  2)  =  J  (y(?),  c,  z)k'(|'.  2)  d?'  +  k'^cc’*’,  2) 


* 


6y 


(27) 


Sr 

where  is  the  functional  partial  derivative,  defined  by 


lim 


R*(Y(C)  -r  ZiC"(C).  C.  z) 

o 


Aiyxcj, 


A. 


„0  ^ 


I_  I  J. 


-00 


R  (Y(C^),  C,  2)  •  k‘(C^,  2)  dC^  (28) 


If  K  is  a  linearized  collision  operator,  we  introduce  a  reflection 

* 

function  R  (u,  V,  w,  u  ,  V  w  ,  z)  such  that 

o  o  o 


-14- 


R*(Y(C„),  Z,  z) 


too  ^ 

1  « 

R  (u,  V,  w,  u^,  v^,  w^,  z) 


X  V(u  ,  V  ,  w  )  du  dv  dw 
o  o  o  o  o  o 


(29) 


and  note  that 


6V 


R  (u,  V,  w,  v^,  w^,  z) 


(30) 


Substituting  Eqa.  (29)  and  (30)  in  Eq.  (27),  we  recover  the  equivalent 
of  Eq.  (9)  : 


s  r*  r**  * 

R  (u,  V,  w,  u  ,  V  ,  w  ,  z)Y(u  ,  V  ,  w  )  du  dv  dw 

02  J  J  J  '  ’  ’  ’  o  o  o  o  o  o  o  o  o 


-  rj 

'j 


.+00  f+<» 


[^r  r  R*^ii  V  u  II '  v'  w' 

J  ^  ""o’  ''o’  u’ 

00  -  00  .  00  O 


X  (u'l  v',  w',  z)  du'  dv'  dw'  +  —  l'^(u,  v,  w,  z)  (31) 
o  o  o  o  o  o  u 


The  initial  condition  is  always  R'(C,  C^,  0)  =  0,  elenming  from  the 

observation  that  there  are  no  collisions  to  reflect  particles  from 

C  to  C  if  T  ■  0. 
o 


-15- 
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